Exact exponents of edge singularities in dynamic correlation functions of ID Bose gas 
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The spectral function and dynamic structure factor of bosons interacting by contact repulsion and 
confined to one dimension exhibit power-law singularities along the dispersion curves of the collective 
modes. We find the corresponding exponents exactly, by relating them to the known Bethe ansatz 
solution of the Lieb-Liniger model. Remarkably, the Luttinger liquid theory predictions for the 
exponents fail even at low energies, once the immediate vicinities of the edges are considered. 
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In a typical experiment with cold atoms, the atomic 
ensemble is kept isolated from environment over certain 
time. This decoupling from the environmental degrees 
of freedom is effective enough to make the interaction 
between the atoms to be the leading cause of the en- 
semble's evolution in time. The isolation of ®^Rb atoms 
in elongated traps recently showed [H the peculiarity of 
the dynamics of bosons confined to one dimension (ID), 
showing almost no relaxation, and consistent with the 
notion of the integrability. 

Integrability allows one to find the spectrum of exci- 
tations of the quantum system by means of a thermody- 
namic Bethe ansatz (TBA), as it was done for the case 
of pointlike interaction between bosons by Lieb A 
more recent progress based on algebraic Bethe ansatz [1] 
ideas helped to develop sophisticated numerical methods 
for studying the dynamic responses of ID Bose gas 0, @| ■ 
However, in spite of more than 40 years of research, the 
analytic calculation of correlation functions from Bethe 
ansatz wave functions still remains a challenge. Some 
progress in understanding of dynamical correlations was 
achieved recently along a different path, which uses the 
effective Hamiltonian description and the physical anal- 
ogy with the Fermi edge singularity [1, 0, H, 0|. The 
analogy helped to show that the dynamic structure fac- 
tor (DSF) and spectral function may exhibit power-law 
singularities along the dispersion curves of the collective 
modes. Very recently Pereira et al. [lo| put forward an 
idea to combine the effective Hamiltonian method devel- 
oped in Refs. [1,0, [i| with the TBA method to study the 
edge singularities of XXZ spin chain. In their approach, 
Pereira et al. used the ideas based on conformal field the- 
ory to relate the parameters of the effective Hamiltonian 
to the finite-size corrections to the energy. 

Here we find an alternative way to directly relate the 
parameters of the effective Hamiltonian to the scatter- 
ing phase shifts found in the Bethe ansatz solution of 
the Lieb-Liniger model. We extend the formalism of 
Refs. 0,0,11 to evaluate the exponents not only of the 
DSF, but of the spectral function as well. In the low en- 
ergy limit we analytically evaluate exact exponents as a 
function of Luttinger parameter K only. Remarkably, the 
Luttinger liquid (LL) theory predictions 11, 13, 3, 14 1 



for the exponents fail even at low energies, once the im- 
mediate vicinities of the edges are considered. Our re- 
sults at low energies help understanding the behavior of 
ID systems beyond the linear spectrum approximation of 
the LL theory. They can be used as a benchmark for nu- 
merical methods which attempt to evaluate many-body 
dynamics of continuous ID models, e.g., using t-DMRG 
algorithms [l5l |. 

DSF describes the probability to excite the ground 
state with momentum and energy transfer (k,Lo), and 
can be measured using Bragg spectroscopy [16| . Spectral 
function describes the tunneling probability for a par- 
ticle (or a hole) with momentum and energy {k,Lu) re- 
spectively. It can be measured using stimulated Raman 
transition combined with additional spin-flips [itI . [i8| . 

In what follows we will be interested in the zero- 
temperature DSF 



S'(fc, tj) = dxdt e*("*-'=^) {p{x, t)p{0, 0)) , 



(1) 



and spectral function A{k,ijj) = — ^ImGffc, uj) signo;, 
where Green's function G{k,u)) is defined by 



G(fc,w) 



dxdie'('^*-'=^) {T (^'(x, t)^'t (0, 0)) ) .(2) 



Here ^'(a;,f) and p{x,t) are boson annihilation and den- 
sity operators, and T denotes time ordering. Energy 
Lu is measured from chemical potential, so A(k,uj) for 
(jj > (w < 0) describes the response of the system to 
an addition of an extra particle (hole). Both A{k, oj) and 
S{k, uj) do not change under transformation k —>■ —k, and 
we will consider them only for fc > 0. Also S{k,uj) ~ 
for a; < 0, so we consider DSF only for oj > 0. 

The exactly solvable Lieb-Liniger model [1, 0| is de- 
fined by 



Hll - 



N g2 



2c 



l<i<fe<Af 



6{z, - zk) - hN. (3) 



Here N is the total number of particles, h is the chemical 
potential, c > is the interaction strength. Hereinafter 
we set h, = 1 and the mass equals 1/2 for brevity. The 
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FIG. 1: (a) Dynamic structure factor (DSF) S{k,Lo) and (b) 
spectral function A{k,u!). Shaded areas indicate the regions 
where they are nonvanishing. Lieb's particle mode ei(fc) and 
hole excitation mode e2{k) are indicated. For spectral func- 
tion A{k,ui) region with uj > {ui < 0) corresponds to the 
particle (hole) part of the spectrum. Notations of jj. indicate 
which exponents of Eqs. (|16|) -(|18 | ) should be used in Eq. Q. 



ground state is fully characterized by the dimensionless 
parameter 7 = c/D, where D = N/L is the density. The 
regime of weak interactions corresponds to 7 <C 1, while 
the strong interactions (Tonks-Girardeau limit) corre- 
spond to 7 ^ 1. An important parameter appea ring in 
the effective hydrodynamic description [ill . hj,hj of the 
model ^ is the LL parameter K = Vf/v, where v is the 
sound velocity and Vf = 27rD is the Fermi velocity of 
noninteracting Fermi gas of density D. Parameter K > 1 
is uniquely defined by 7, with K « 7r7~^/^ for 7^1 and 
if « 1 4/7 for 7 > 1 0, [il. Energies of excitations 
can be also obtained, and give rise to Lieb's particle (ei) 
and hole (£2) excitations, as shown in Fig|T] The mini- 
mal excitation energy (measured from the ground state) 
of a state with total momentum k equals £2 (A:), and both 
A{k^uj) and S{k,uj) vanish identically for \u!\ < e2{k). 

At zero teniperature, S{k,uj) and A{k,uj) have power- 
law behavior [a, 0, S| near ±£1(2) (fc), 



S{k, oj), A{k, iv) ~ const -I- 



UJ ± £i(2)(fc) 



(4) 



with notations of /i shown in Fig. [TJ We provide exact 
results for exponents at ±£i(fc) with fc > 0, and for expo- 
nents at ±£2(fc) with < fc < 2TrD, see Eqs. ^ - ([^ . 

Let us briefly review the solution of the Licb-Linigcr 
model to set the notations, and we will mostly follow 
the conventions of Ref . 3] . Ground state quasimomenta 
VjT^ < i < N, are given by the solutions of Bethe equa- 
tions 



JV 



Li/j + 0{uj — Uj^) — 2TTnj, 



(5) 



*:=1 



where 9{x) = 2 arctan ^ is the two-particle phase shift 
and quantum numbers are Uj = J — 1 — {N — l)/2. In 



the thermodynamic limit, this system gives rise to the 
integral equation 



K{v, ii)p{fi)dfi — 1. 



(6) 



Here p{v) = lim 1/ (L{vk+i — i^k)) is the density of roots, 
K{i^, p) = , and ±q is the highest (lowest) filled 

quasimomentum; q is defined as a function of density by 
the normalization condition D = J^^ p{i')di'. Particlelike 
excitations with fc > can be constructed by adding an 
extra quasimomentum X > q, while holelike excitations 
are obtained by removing a quasimomentum |A| < 5 (par- 
ticlelike states with fc < correspond to A < — g). Since 
all quasimomenta Vj are coupled to each other by Eqs. 
([5]), this will shift all of them. A convenient way to take 
this change into account is to introduce a shift function 



(7) 



where Vj are new solutions and upper (lower) sign cor- 
responds to extra particle (hole) . In the thermodynamic 
limit Fb{v\\) satisfies an integral equation [IBl 

FB{y\\) ^ I' K{v,p)FB{p\X)dp = ^±fc^.(8) 

Shift function Fb(i/|A) can be used 's'l to calculate exact 
energies £1(2) ^ a function of momentum fc(A). The latter 
can be also written Q as 



fc(A) = ± (^A - ttD + y""^ 



Q{\ - v)p{v)dv 



(9) 



Here upper (lower) sign corresponds to particle (hole) 
excitation with X > q (|A| < q), and k{q) = 0, k{—q + 
0) = 27rL>. 

As will be shown below, FB{±q\X) play a crucial role 
in the calculation of the edge singularities, so we will 
investigate it in more detail. One can analytically derive 
the limiting behavior 

FB{±q\q) = 1/2 ± (1/2 - 1/{2Vk)) , (10) 

and FB{±q\ — q) ~ \fK — FB{^q\q)\ in addition, 

Fi3(±(j|A) « c\/X/(7rfc) < 1 for q,c<A«fc. (11) 

Eq. (dni) can be derived ^ from Ref. [lH. Eq. HID 
follows from expansion of right hand side of Eq. ([5]) com- 
bined with p{±q) = \/~K/2tt, see e.g. Eqs. (I.9.20-I.9.22) 
of Ref. i. 

We can calculate the exponents using the method of 
Refs. 0,0, S, For example, exponents at £i(fc) > 
are evaluated using effective Hamiltonian H = Hq + Hd + 
Hint, where 



H^ = ^jdx (^K{Vef + ;^(V0)2^ , (12 
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Hd= I dx d\x){e,{k) ~ ^^^^)d{x), (13) 



H^nt ^ dx[ Vr\/ 



27r 



VlV- 



2tt 



dHx)d{x).{U) 



Here creates an extra particle with momentum near 
k, and we use the notations of Ref. [iJIj such that 
-^Bix.t) ~ e*^(^'*\ [(j){x),V9{x')] = iirSix - x'). The ef- 
fective Hamiltonian, Eqs. P^ - p^ . discriminates states 
created by d'l' from low-lying states of the rest of the sys- 
tem which are described as LL. This separation is pos- 
sible only in the investigation of response functions very 
close to the Lieb's modes, such as line ei{k) in Fig. [2l 
and only due to the finite curvature of the underlying 
boson spectrum. Under these circumstances, excitations 
close to the edge can be distinguished from the excita- 
tions of the LL. Response functions evaluated from Eqs. 
(I12p - (|14l) are valid in a narrow region of width vanishing 
as ~ fc^ at fc ^ 0; see discussion after Eq. (|27|) . 
Singular parts of DSP and spectral function are 

Sik,uj)^ J dxdte'^'{d^ix,t)^d^iO,0))„„+H,+H,„„ 

A{k,uj)^ [ dxdte'^'{dix,t)dHO,0))Ho+H,+H„,,- 



The crucial step in our approach is the identification of 
the unitary transformation which removes an interaction 
term Hint from the Bethe ansatz solution. This can be 
done by noticing that Hq becomes a sum of noninteract- 
ing modes after transformation = \/K(f),6 = 
If one refermionizes fields 0, 9, one obtains a noninter- 
acting Luttingcr model with two branches. Fermionic 
excitations of this model with momenta kj > corre- 
spond to low-energy particlelike excitations of the Bethe 
ansatz with quasimomcnta Vj > q, and kj oc Vj — q. 
The state where a d particle is present corresponds to a 
state with occupied quasimomentum A. For noninteract- 
ing fermions, phase shift on the particle d can be written 
as 6{k) — 27r(fcj — fcj)/(/cj+i — kj), where kj is the mo- 
mentum of the new eigenstate in the presence of particle 
d. The same quantity in the Bethe ansatz solution can be 
calculated using Eq. ([7]) , and by noticing that the shift of 
quasimomenta by I'j+i — Vj corresponds to a phase shift 



27r. Then interaction term Hint can be removed [23l. |24| 
by unitary transformation U^{Ho + Hd + Hint)U, where 

S±=2TTFB{±q,X). (15) 



A standard calculation [IJ, [23|, |24|, [251 then leads to 
1 / 1 



Ml,2 = 1 



K 



2tt 



— 1 1 f ^+-^- 



2 V 27r 

1 / 6++S^ 

2 I 2n 



,(16) 
(17) 



Exponents /i+ and /i_ are more complicated. For 
example, ^+ corresponds to a state with one addi- 
tional particle, total momentum fc, and the smallest 
possible energy £2{k). Such state is given by two ex- 
tra particles at the right quasifermi surface and a hole. 
To evaluate one needs to calculate the correlator 
(d'l'vI/\I/(x,t)vI/tv[/t(f(o, 0)), where d creates a hole. One 
obtains 



1 



K 



27r 



2ti 



(18) 



In Eqs. (|T6|) - p8p one should choose A > (7 for /ii, 
and |A| < q for /X2,/i+,7I^. 

Using Eqs. (nnD-ini), one can derive analytically the 

behavior of exponents near k = 0, 27rZ) and at fc — > 00. 
DSF exponents are given by 

/ii,2(0) = 0, (19) 
Ijl2{2'kD ^Q) = 2^/K -2K <Q, (20) 
/ii(oo) = 1 - l/(2ii:) > 0. (21) 

Exponents for the spectral function at low energies equal 



71^(0) = 1/s/K - 1/{2K) > 0, (22) 
M±(0) = -1/Vk -1/{2K) < 0, (23) 
7iT(27rL>- 0) = -2X + 2\/X + l/VX- l/(2i^), (24) 
^i+{2nD - 0) = -2K + 2Vk - 1/Vk - 1/{2K), (25) 



and for large momenta one has 

/i_(oo) = 1 - 2/K. 



(26) 



For k CO Eq. PT)) gives /i+ ^ 1, which is a certain way 
to parameterize an expected result A(k,u!) ^ S{uj — fc^), 
since S{uj - k^) ~ limjr;:^! / dte'^'^-''^'>yt^'~ . 

Numerical solutions show that all exponents are mono- 
tonic functions of k, lying between limiting values cited 
above. Exponents JT^ and /i_ can change sign as a func- 
tion of k. Exponent changes its sign from positive at 
/c = to negative at A; = 2TrD for sufficiently weak inter- 
actions, a,t K > 1/4-1- V5/4-l-l/2y^l/2-|-\/5/2« 
1.445, and can also change its sign from negative at 
fc = to positive at /c = 00 for X > 2. Our results 
are in full agreement with previously published limiting 
cases 0, and agree qualitatively with the numerical 
calculations of DSF ^ and spectral function . 

One of the main achievements of the LL approxima- 
tion [11], [13, [13, 14, 25 1 is the calculation of exponents 
in dynamic correlation functions in the limit k —>■ 2TTDn, 
where n is any integer. We note however that exact ex- 
ponents ([20|) - (l25p in the immediate vicinity of the edges 



are different. They show markedly non-LL behavior, e.g. 
they can depend on ^/K, while LL exponents depend only 
on K and 1 /K. To illustrate the failure of LL exponents 
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Luttinger 
liquid 




FIG. 2: Spectral function A{k, to) for uj > near k — 0. Exact 
exponents TTj: = - 1/{2K), fi^= -XjyfK -Xj {2K) are 

valid only in the immediate vicinities of ei and £2 marked by 
dashed lines. Luttinger liquid behavior given by Eqs. (|27|) 
is valid only sufficiently far away from e = life, for uj — vk"^ 
ei - £2. 



in the immediate vicinity of the edge we consider A(k, lS) 
for a; > near A: = in more detail; see Fig. [2l LL 
theory assumes linear spectrum e(fc) ss vk and predicts 
H, [iJI at fc <C u 



0, if Lo < vk, 
{oj - t>A:)i/(4A')-i if < w - w/c < w, (27) 

^1/{2K)^2 if ^fc <^ ^_ 



Difference between the second and third lines here arises 
from the fact that only one branch contributes to the 
exponent for < to — vk ^ uj, while both left- and right- 
movers contribute at vk <^ lo. One should note however, 
that in condition < — f fc <C cj, is an artifact of the 
linear spectrum approximation, which breaks down at 
energy scales on the order of ei{k) — e2(fc) ~ fc^. It is pre- 
cisely within this region where exact results (|22p . ((23| are 
applicable. Each of these exponents is valid in the vicini- 
ties of £1(2) (fc) much narrower than ei{k) — £2{k) ~ fc^. 
These exponents describe the response of the system at 
fc — > beyond linear hydrodynamic approximation. Nev- 
ertheless, for the Lieb-Liniger model they turn out to de- 
pend only on the LL parameter K, since phase shifts at 
quasifermi surface are also related to K] see Eq. ([T0| . 
Far away from k — > 27rDn, regions of validity of expo- 
nents (fT7)) . p^ widen, and the singularities constitute the 
main features of the spectral function. 

To summarize, we have found exact exponents char- 
acterizing the singularities of the dynamic structure fac- 
tor and spectral function along the dispersion curves of 
the collective modes of a ID Bose gas. The found expo- 
nents are related to the known Bcthe ansatz solution of 
the Lieb-Liniger model. Remarkably, the Luttinger liq- 
uid theory predictions for the exponents fail even at low 
energies, once the immediate vicinities of the edges are 
considered. 
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Note added. - After the essential par t of this work 
has been completed, preprint to Ref. [26| has appeared, 
where equations similar to our Eqs. ([8]) and (|16p were 
also reported for DSF of ID fermions on a lattice. 



[1 
[2: 

[s: 

[4] 
[5] 

[6: 

[7 

[s: 
[9: 
[10: 

[11 

[12 

[13 
[14 

[15 
[16 

[17 
[18 

[19 
[20 



T. Kinoshita, T. Wenger, and D. S. Weiss, Nature 440, 
900 (2006). 

E.H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963); 

E. H. Lieb, ihid. 130, 1616 (1963). 

V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin, Quan- 
tum Inverse Scattering Method and Correlation Func- 
tions (Cambridge University Press, Cambridge, 1993). 
J.-S. Caux and P. Calabrese, Phys. Rev. A 74, 
031605(R) (2006). 

J.-S. Caux, P. Calabrese, and N.A. Slavnov, J. Stat. 
Mech. (2007) P01008. This work considered A{k,uj < 0). 
M. Pustilnik, M. Khodas, A. Kamenev, and L.I. Glaz- 
man, Phys. Rev. Lett. 96, 196405 (2006). 
M. Khodas, M. Pustilnik, A. Kamenev, and L.I. Glaz- 
man, Phys. Rev. B 76, 155402 (2007). 
M. Khodas, M. Pustilnik, A. Kamenev, and L.I. Glaz- 
man, Phys. Rev. Lett. 99, 110405 (2007). 
M. Khodas, A. Kamenev, and L.I. Glazman, 
arXiv:0710.2910vl. 

R.G. Pereira, S.R. White, and I. Affleck, Phys. Rev. Lett. 
100, 027206 (2008). 

K.B. Efetov and A.I. Larkin, Sov. Phys. JETP 42, 390 
(1975) . 

F. D.M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 
M.A. Cazalilla, J. Phys. B 37, SI (2004). 
T. Giamarchi, Quantum Physics in One Dimension 
(Clarendon, Oxford, 2004). 

U. SchoUwock, Rev. Mod. Phys. 77, 259 (2005). 
D.M. Stamper-Kurn et al, Phys. Rev. Lett. 83, 2876 
(1999). 

L.-M. Duan, Phys. Rev. Lett. 96, 103201 (2006). 
T.-L. Dao, A. Georges, J. Dalibard, C. Salomon, and 
I. Carusotto, Phys. Rev. Lett. 98, 240402 (2007). 
A. A. Abrikosov, L.P. Gorkov, and I.E. Dzyaloshinski, 
Methods of Quantum Field Theory in Statistical Physics 
(Dover, New York, 1963). 

Shift function F{u\X) defined by Eq. (1.4.25) of Ref. 0] 
is related to ours as F{i'\X) — Fb{i^\X) — -npiv). Ref. 
[3| assumes a change of periodic boundary conditions to 
antiperiodic ones with the change of particle number (see 
p. 23), and phase shifts calculated there correspond to the 
fermionic Cheon-Shigehara model [IJl dual to the bosonic 
Lieb-Liniger model. 
[21] T. Cheon and T. Shigehara, Phys. Rev. Lett. 82, 2536 
(1999). 

[22] V. Korepin and N. Slavnov, Eur. Phys. J. B 5, 555 (1998). 
[23] K.D. Schotte and U. Schotte, Phys. Rev. 182, 479 (1969). 
[24] L. Balents, Phys. Rev. B 61, 4429 (2000). 
[25] A. Gogolin, A. Nersesyan, and A. Tsvelik, Bosonization 

and Strongly Correlated Systems (Cambridge University 

Press, Cambridge, 1998). 
[26] V.V. Cheianov and M. Pustilnik, Phys. Rev. Lett. 100, 

126403 (2008). 



